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Granular material on an inclined plane will flow like a fluid if the angle 9 the plane makes with the 
horizontal is large enough. We employ a modification of a hydrodynamic model introduced previ- 
ously to describe Couette flow experiments to describe chute flow down a plane. In this geometry, 
our model predicts a jammed-to-flowing transition as 9 is increased even though it does not include 
solid friction, which might seem necessary to stabilize a state without flow. The transition is driven 
by coupling between mean and fluctuating velocity. In agreement with experiments and simulations, 
it predicts flow for layers with a thickness H larger than a critical value H stop {9) and absence of 
flow for H < // s top (6*). 



In recent years, considerable experimental and theoret- 
ical effort has been devoted to the study of flow properties 
of granular media, yet no fully satisfying phenomenology 
of granular rheology has emerged jij. The major diffi- 
culty in describing granular flows, particularly in dense 
media, is that both kinetic phenomena associated with 
particle motion and more static phenomena like solid fric- 
tion and force chains are expected to be important. Fol- 
lowing early work of Jenkins and Savage and of Haft , 
various kinetic-hydrodynamic threories based on coupled 
equations for the transport of mcmentum and kinetic en- 
ergy have been constructed. They assume instantaneous 
binary collisions and neglect any solid-like friction be- 
tween grains. Various authors have proposed including 
shear independent terms in the stress tensor to capture 
empirically the effects of intergrain friction Nei- 
ther of these types of theories incorporates the effects 
of strongly inhomogenoeus force networks in static and 
slowly driven dense granular materials |||]] 

On the basis of experimental results on sheared granu- 
lar material in Couette geometry, we recently proposed a 
modified hydrodynamic description Q], whose aim is to 
include high-density effects in equations describing en- 
ergy and momentum transport. Specifically, in analogy 
with the behavior of supercooled liquids close to the glass 
transition, we assumed an "anomalous" divergence of the 
viscosity close to random close packing. Even though our 
approach still does not include enduring contacts and 
solid friction between grains, it is able to reproduce most 
of the experimental results in the Couette Geometry. 

Granular flow down inclined slopes is of considerable 
practical importance in geophysical phenomena such as 
rock and snow avalanches, pyroclastic flows, etc.. Sys- 
tematic experiments m and molecular simulations fic[ | 
on a layer of beads on an inclined plane as a function 



of layer thickness H and angle of the plane relative 
to horizontal yield a curve H stop (9) separating jammed, 
stationary behavior at small H from flowing behavior at 
large H . Moreover, precise and systematic measurements 
of the mean-bead velocity yield a maximum velocity that 
scales as H 3 / 2 . 

In this paper, we analyze chute flow down inclined 
planes using the frictionless kinetic-hydrodynamic the- 
ory developed for the study of Couette flows 0. Our 
analysis produces results in agreement with experiment 
and simulation in this geometry and provides further ev- 
idence that our phenomenological kinetic-hydrodynamic 
theory provides a robust and general description of flow- 
ing granular media. In particular, we find an h stop curve, 
or equivalently a critical angle separating stationary from 
flowing behavior. Normally, the existence of such a criti- 
cal angle 6 C is associated with solid friction: a solid block 
on an inclined plane is stationary for 9 < 9 C and slides for 
9 > 9 C . In our analysis, the existence of C is associated 
with a purely dynamical jammed-to-flowing transition in 
which solid friction is totally absent. The origin of this 
transition is the coupling between velocity fluctuations 
and mean flow - an important feature of our previous 
study. 

Geometry of the problem. We consider the geome- 
try of the experiments in Refs. Q : a layer of granular ma- 
terial, composed of spheres of diameter d, of thickness H 
flows on an inclined plane making an angle 9 with the hor- 
izontal. We choose a coordinate system with the x-axis 
parallel to the plane, the z-axis perpendicular to it, and 
origin at the bottom of the flowing layer. In a stationary 
state, conservation of momentum in directions parallel 
and perpendicular to the plane leads to the equations: 
d z a xz — p(z)gsm9 and d z a zz — p(z)g cos9, where p is 
the mass density and a zz and a xz are, respectively, the 
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diagonal and off diagonal components of the stress ten- 
sor. Combining these two model-independent equations 
yields the simple relation: a xz — ta,n6<j zz + C, where C is 
a constant independent of z. In our hydrodynamic model, 
the off diagonal component of the stress tensor is viscous 
like: o xz — r)(p, T)i, where 7 = dV x /dz is the shear rate, 
V x the mean velocity, and r\ is the viscosity, a function of 
the density and granular temperature field T. As usual 
the latter is defined in terms of the RMS part of the ve- 
locity field T(z) = ^m(v 2 (z)}. The temperature profile 
is determined by the balance between viscous heating, 
heat flow, and dissipation through inelastic collisions |Q : 
£\(p,T)JLT + a xz j-e(p,T)T = 0, where X(p,T) is the 
thermal conductivity and e(p, T) is the rate of kinetic en- 
ergy loss. These hydrodynamic equations are closed by 
the equation of state P = f(p)T. In contrast to classi- 
cal thermalized fluids, there arc strong temperature and 
density variations over the system, and the density and 
temperature dependence of transport coefficients play an 
important role in determining flow properties. These de- 
pendences are usually obtained using kinetic theory ||] 
with the Enskog approximation for the collision kernel, 
which accounts for excluded volume but neglects any cor- 
relation between collisions. By construction this theory 
yields transport coefficients that are roughly proportional 
to the collision frequency, i.e., to the pair-correlation 
function at contact, which in the high-density limit is 
expected to diverge like (p c — p)^ 1 close to the max- 
imum density, p Cl allowed by excluded volume effects 
jnj. The assumption of negligible correlation between 
collisions should fail at high density, and in our mod- 
ified description of Ref. ]7j], we proposed replacing the 
Enskog expression for the viscosity with one with a spe- 
cific "anomalous" algebraic divergence, r\ ~ (p c — p)^ 13 
fl3f close to p c , in full analogy with the behavior of su- 
percooled liquids close to the glass transition j^J. This 
picture is supported by the experimental observation in 
sheared granular material of a specific scaling law relat- 
ing the local shear rate 7 to the granular temperature T 
IP- 

In the Couette geomety, the density is high through- 
out the cell. In contrast, in the current chute geome- 
try, small densities could in principle be reached at the 
top surface, and we need to include the low-density limit 
in our expressions for transport coefficients. Since our 
aim is to provide a sketch of the jammed-to-flowing, we 
shall not use the full Enskog expression for transport co- 
efficient in the small and intermediate density regions, 
as is done by various authors Instead, we will use 
simple interpolated expressions that are consistent with 
both the high-, intermediate-, and small-density limit. 
Our specific choices are X(p,T) = [A /(to 1 ' 2 d 2 ))(l — 
p/pc)- 1 ^ 2 , e(p,T) = [ev/( m ^ 2 d)]p(l- p/pc)- 1 ^/ 2 , 
rj{p,T) — rjolm 1 / 2 / d 2 ]g{p)T 1 / 2 , where m the mass of 
one particle, d is its diameter, and Ao, £q and 770 are 
dimensionless constants. The function g{p) is chosen 



so as to interpolate between the very low-density limit, 
where the viscosity goes to a constant, to its high- 
density "anomalously" diverging behavior close to p c , 
passing through its Enskog behavior for intermediate 
(but already high) densities. A sensible choice is sim- 
ply g(p) = {a{pf- 1 + a^ 1 )a(p)^/(l + c^" 1 ), where 
a(p) = 1 — pj p c and (3 is the anomalous exponent dis- 
cussed above. The parameter ao fixes the density where 
the crossover to the anomalous regime occurs. The equa- 
tion of state is written as P/T = p/(l — p/p c ), which 
does indeed reduce to the correct limits in the low- and 
high-density limits O]. We emphasize that our results 
are not very sensitive to the specifics forms for g and 
the equation of state. Our expressions were chosen to 
minimize the number of parameters in the model. 

The previous equations can be made dimensionless us- 
ing a reduced velocity v — V/ \/{gd) , temperature t = 
T/{mgd), stress a = o xz /(p c gd), pressure p = P/{p c gd), 
and distance z — z/d (we shall omit the tilde over z in the 
following). Using the previously defined expressions for 
the transport coefficients, we can now rewrite the coupled 
equations for stress and temperature as: 

£P = - C0sd ^r t ( la ) 
d z d z t) + a g{ _£ )tl/2 - bpt^ 2 = 0, (fb) 

where p/(p + t) = p/p c is the reduced density, a = 
(p c d 3 ) 2 /(X0V0) an d b = e p c d 3 /Xo are dimensionless con- 
stants, a is related to p according to a = tan Op + c, with 
c an as yet unspecified constant. The parameters ap- 
pearing in Eq. (Q) are strongly constrained by the com- 
plementary experiments performed in the Couette cell 
|^,^5|. In this case, the velocity and temperature profiles 
are localized close to the moving boundary over a depth 
of a few particles diameter. This behavior is predicted 
by the hydrodynamic model with a characteristic decay 
lentgh for the RMS velocity given by S/d = (2/fe) 1 / 2 @. 
Experimentally S is of the order of 4-5 particle diameters, 
yielding b ~ 0.1 (we shall use in the following b — 0.111). 
The exponent [3 was determined experimentally in Ref. 
|Q to be 1.75, while experiments over a larger range of 
data by Mueth yielded (3 = 1.5. The parameter ao is 
chosen such that the cross over to the anomalous scal- 
ing occurs roughly 10% below random close packing. In 
pratice ao = 0.05 is adequate. Eventually, a is expected 
to be of order unity. 

The hydrodynamic equations (jl]) have to be supple- 
mented by boundary conditions for the velocity and tem- 
perature fields at the bottom and top surfaces. As to the 
conditions on the granular temperature, we showed in our 
previous study of the Couette cell that a vanishing "heat 
flux" condition at the boundaries yielded results in agree- 
ment with the experimental ones. In the present study, 
we shall thus assume d z T = at the boundaries. In the 
recent simulations of Silbert et al. |10[| , the temperature 
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profile has a vanishing derivative close to the bottom and 
top boundaries. The previous assumption thus amounts 
to slightly shifting the hydrodynamic boundaries of the 
system, a behavior which is encountered in classical flu- 
ids The boundary condition for the mean velocity 
field are a priori more obvious: at the wall surface, the 
velocity is assumed to vanish, while at the top surface a 
stress-free boundary condition will be assumed, a = 0, 
yielding d z v = at z = H. Since, as mentioned above, 
the effective hydrodynamic boundary z = H is located 
slightly inside the material, the pressure is not expected 
to vanish there, and it can be written as po = z o cos 9 
(in reduced units) with zq a molecular distance of the or- 
der of the diameter (unity in reduced units). To be fully 
consistent with the stress free boundary for the velocity, 
we thus write a = tanf?(p — po) (i.e., fix the constant 
c = — tanOpo). As emphasized above, this is a molecu- 
lar effect that disappears in the large- H limit. Moreover 
this choice does not affect the existence of the jammed- 
to-flowing transition to be discussed now. 

The previous system of equations (|l|) has been solved 
using a standard Runge Kutta algorithm, using a 
shooting-like procedure: for a given temperature T(H) 
at the top boundary, the derivative of T at the bottom, 
d z T(0) is computed. The solutions we seek thus corre- 
spond to zeros in the d z T(0) versus T(H) curve. A tran- 
sition between two kinds of behaviour is then found: for a 
set {H, 9} such that H is below a critical H stop (9) curve, 
only the T = solution exists. This corresponds to a 
"jammed" state, which stays at rest. Above this critical 
value, nonzero solutions for the temperature, which cor- 
respond to a flowing regime do exist. In Fig. [j], we plot 
the boundary between the two regimes, denoted H stop (9) 
as in Ref. || . This curve is in qualitative agreement with 
the experimental results of Refs pi . 




FIG. 1. Phase diagram for the jammed-to-flowing transi- 
tion in a gravity driven granular material as obtained from 
the resolution of the hydrodynamic equations (|l|). The pa- 
rameters used are a — 1.3, b — 0.111, Qo = 0.05, zo = 0.5. 

This shows that the hydrodynamic model is able to 
generate a finite critical angle below which no flow oc- 
curs. Within the model, the origin of this behavior is 



simply the balance between the "viscous heating" , which 
generates fluctuations, and dissipation through collisions, 
which tends to inhibit the flow. The transition is thus 
purely dynamical: a jammed state occurs for low an- 
gles/thicknesses because fluctuations are insufficient to 
allow flow. 

Above H s to P (9), more than one solution can be found 
for some sets of parameters {H, 9}. However a linear 
stability analysis of the hydrodynamic equation for the 
temperature shows that only a single nonzero solution 
is dynamically stable. On the other hand the T = 
solution is always linearly stable but becomes unstable 
against larger "kicks" in the temperature. It is impor- 
tant to note that there is a solution to the hydrodynamic 
equations for any set of H and 9 above the boundary 
H s top(9), as is found experimentally. This is in contrast 
to other approaches where for a given 9 a single flowing 
thickness H is selected |J. In the flowing regime, we ob- 
tain temperature and density profiles, which are plotted 
in Fig. | for H = 40. 
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FIG. 2. Top : Normalized temperature profiles for H = 40, 
6 = 20.55, 20.90, 21.50, 22.50. Bottom : Normalized density 
profiles for H = 40, 9 = 20.55, 20.90, 22.50, 24.50, 26.50 from 
top to bottom, (a = 1.3, b = 0.111, «o = 0.05, zo = 0.5) 

Both results for the density and temperature are in 
qualitative agreement with the simulations of Refs. [p~0|| . 
In particular the density profile is found to be almost 
constant in the middle of the sample, with a large drop at 
the top surface, as found numerically. In the simulations 
however, the plateau in the density is flatter than ours. 

Velocity profiles are obtained by integrating the shear 
rate j = dv x /dz, which in our model is linearly related to 
the non diagonal stress cr xz (z) — r]{p{z),T{z)) A f{z). We 
find that the z dependence of the profiles is in very good 
agreement with the predictions of a simple Bagnold scal- 
ing a = A Bag 7 2 , vW/vmax = 1 - (1 - z/H) 3 / 2 . Within 
the Bagnold assumption, v max = (2/3) AeagVsin 9H 3 / 2 . 
This agreement simply indicates that once the system 
flows, Bagnold-like scaling is a good approximation to the 
flow properties, even though it cannot predict a jammed- 
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to-flowing transition. The same agreement with Bagnold 
scaling was observed in the simulations of Ref. . The 
_ff 3 / 2 dependence of the maximum velocity t> ma x is also 
recovered in our model in agreement with experiments 
and simulations. Our results are also consistent with a 
V sin 9 dependence of t/ max with v4e ag ~ 0.2. This is in 
fair agreement with the simulations of Silbert et al. [jl0| , 
who find a slight dependence of ^Bag on 9. On the other 
hand, the experiments of Pouliqucn exhibit a stronger 
dependence on 9, with v max scaling like [H s t op (9)]^ 1 

Discussion. Reference M introduced a modified 
kinetic-hydrodynamic theory for granular flow that pre- 
dicted temperature and velocity profiles in dense systems 
in Couette geometry in agreement with experiments. In 
this paper, we calculate the corresponding profiles for 
flow down an inclined plane using a generalization of this 
theory to include gravity and low- as well as intermediate- 
and high-density regimes. Our results agree with those 
obtained in experiments || and simulations fl(i[| . Thus, 
the kinetic-hhydrodynamic theory of granular flow pro- 
vides a good description of flow in two qualitatively dif- 
ferent geometries. Though it is premature to generalize 
from two to all geometries, it would, nevertheless, ap- 
pear that this description might serve as an all purpose 
phenomenological theory for granular flow. In the case 
of flow down a plane, it is able, contrary to intuition, to 
predict a jammed-to-flowing transition and a critical an- 
gle below which no flow takes place, even though it does 
not include solid friction between grains, which might be 
viewed as necessary to stabilize a state without flow. 

Our results raise a number of basic questions about 
the mechanisms underlying granular flows. Two impor- 
tant features of granular materials are (i) the possibility 
of solid friction between grains and (ii) the absence of ve- 
locity fluctuations (i.e., temperature) when the system is 
at rest. Under flow, the effects of these features become 
intertwined, making it difficult to provide first principles 
modeling of the rheology of granular materials. The fact 
that our frictionless hydrodynamic model produces suc- 
cessful predictions for two different geometries suggests 
that solid friction is not necessarily a dominant feature 
(though it may not be negligible). On the other hand, 
our theory does highlight the importance of coupling be- 
tween velocity fluctuations and mean flow. Intuitively, 
one might expect this to be a generic feature of rheology 
of "T = 0" glassy systems: in order to flow, a system 
must create its own fluctuations, which in turn couple 
to flow. It is interesting to note that various descrip- 
tions based on similar ideas have recently emerged in at- 
tempts to provide a phenomenological description of the 
rheology of complex systems such as gels, paste 17 , and 
foam |L8|. The hydrodynamic model successfully cap- 
tures coupling between mean and RMS velocity fields, 
and it should serve as a basis for developing improved 
models that could for example include contact forces. At 
high density, enduring contacts occur, but they must not 



be confused with the existence of solid friction between 
grains. Indeed, enduring contacts do not invalidate a hy- 
drodynamic description; they can be taken into account 
in transport equations (as done for example in simple liq- 
uids with continuous interactions (l6|). In this case, only 
the T 1 / 2 scaling of transport coefficients, which stems 
from the assumption of binary collisions, will be mod- 
ified, and a more complex temperature dependence of 
transport coefficients is expected. 
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